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ABSTRACT 

- This  thesis  addresses  the  design  and  analysis  of  discrete 
lateral  autopilots  for  application  to  BTT  missiles. 

The  first  part  review^4. the  classical  design  and  analysis  of 
the  continuous  uncoupled  yaw  and  roll  channels,  ^s  developed  in 
[Ref.  6].  Then,  applying  analog-to-digital  conversion,  the 
corresponding  discrete  autopilots  were  designed  and  analyzed  in 
terms  of  their  transient  responses. 

The  second  part ;,dtilize^  modern  control  design  techniques  for 
the  single-input  discrete  lateral  autopilots.  At  first,  assuming 
availability  of  all  states  for  feedback  purposes,  a  discrete 
state-feedback  autopilot ^ was  obtained.  Next,  since  the  state 
vector  is  not  always  available  to  direct  measurement ,  ,-ap  estimator 

I 

was  introduced  to  implement  control.  The  state-feedback  and 

estimator  designs ,  were  analyzed  for  both  lateral  channels  and 

found  to  have  satisfactory  time  responses. 

Finally,  coupling  the  discrete  pitch  and  roll  channel 

'•  ■  c 

autopilots,  a  state-feedback  and  estimator , werevdesigned  and 

i 

found  to  be  robust,  / 
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I.  INTRODUCTION 


Modern  tactical  missiles  require  increased  stand-off  ranges 
and  need  to  meet  threats  from  highly  maneuverable  air  targets. 

The  high  maneuverability  of  air  targets  has  directed  the  use  of 
defence  missiles  capable  to  develop  higher  lift  accelerations  and 
more  complex  control  laws.  In  order  to  accomplish  the 
requirement  of  large  stand-off  ranges,  propulsion  systems  using 
air-breathing  engines  have  been  studied  and  developed  in  recent 
years.  The  advent  of  air-breathing  engines  has  naturally  led  to 
the  consideration  of  BTT  missiles  in  order  to  minimize  the  inlet 
angle-of-attack. 

The  necessity  of  more  complex  control  laws  has  introduced  the 
application  of  modern  control  and  estimation  theory,  since  more 
complicated  information  of  the  missiles  states  are  needed.  BTT 
controlled  missiles  are  generally  characterized  by  increased 
maneuverability  and  considerable  drag  reduction  over  conventional 
cruciform,  roll  stabilized  STT  controlled  missiles.  Certain 
limitations  in  technology  [Ref.  1]  have  delayed  the  development 
of  BTT  control  systems,  and  consequently  any  progress  in  the  area 
of  BTT  autopilots. 


Major  technological  improvements  during  the  last  decade,  as 
the  availability  of  advanced  digital  computers,  reopened  the 
issue  and  made  BTT  control  feasible  in  spite  the  added  complexity 
of  control  laws  for  the  autopilots.  In  addition,  certain  types 


of  ramjet  engines  [Ref.  2],  which  are  candidate  propulsion 
systems  for  modern  tactical  mission  requirements  of  range  and 
high  altitude  [Ref. 3],  have  presented  a  need  for  a  missile  control 
technique  to  maintain  effective  inlet  flow.  This  was  the  main 
reason  for  given  further  impetus  to  the  investigation  and 
development  of  BIT  control. 

Despite  the  fact  that  BTT  steering  may  provide  improved 
performance  for  a  missile  system,  there  are  still  unanswered 
questions  concerning  stability  during  homing  phase,  guidance 
performance,  autopilot  guidance  logic  and  subsystem  requirements. 
All  these  questions  have  to  be  investigated  and  properly  answered 
in  order  for  BTT  steering  to  be  considered  as  a  viable  control 
method  for  high  performance  missiles. 

During  the  past  decade  many  missile  programs  [Ref.  3j  were 
initiated  to  improve  the  capability  of  steering  tactical  missile 
via  BTT  control  with  results  that  have  greatly  advanced  the 
understanding  of  the  various  missile  subsystems.  In  the 
autopilot  area  many  different  types  have  been  designed  and 
developed.  All  of  them  force  the  missile  to  roll  or  bank,  so  that 
the  steering  maneuver  occurs  with  the  missile  axis  oriented  in  a 
specific  or  preferred  direction  with  respect  to  the  incoming 
airstream.  This  class  of  autopilots  is  usually  known  as  POC 
autopilots . 

The  main  criterion  for  the  selection  of  a  particular  type  of 
autopilot  is  based  upon  the  guidance,  airframe  and  propulsion 
system  requirements.  Generally,  missiles  with  either  one  or  two 
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planes  of  symmetry  use  a  POC  autopilot  which  forces  the  missile 
to  bank  in  order  to  turn  as  an  aircraft.  If  this  motion  is 
coordinated,  then  the  autopilot  is  referred  to  as  a  CBTT 
autopilot . 

In  the  guidance  area,  radome  aberation  effects  for  frequency 
guidance  are  of  major  concern  [Ref.  3]  and  are  being  investigated 
in  great  extent.  Also,  the  interaction  between  BIT  control, 
antenna  stabilization  and  sensor  orientation  are  some  of  the 
additional  concerns  that  have  to  be  properly  addressed.  However, 
simplified  studies  [Ref.  4],  which  neglect  radome  effects  and 
assume  that  the  missile  motion  is  entirely  coordinated,  have 
proven  that  BIT  control  can  provide  acceptable  performance  with 
roll  rates  that  are  not  excessive  for  autopilot  design.  These 
studies  were  made  for  a  medium  range  area  and  long  suppression 
mission,  and  considered  both  high  lift  (i.e.,  planar)  and 
moderate  lift  (i.e.,  cruciform)  airframe  configurations. 

In  order  to  take  full  advantage  of  CBTT  control,  planar 
airframes  have  been  designed  to  increase  the  lifting  capability 
in  one  direction  without  the  weight  and  drag  penalty  associated 
with  orthogonal  lifting  surfaces  [Ref.  5].  These  airframes  have 
aerodynamic  properties  characterized  by  increased  potential  to 
enhance  CBTT  control. 

The  present  thesis  addresses  the  design  and  analysis  of 
discrete  lateral  autopilots  for  application  to  CBTT  missiles. 


The  first  part  reviewed  the  design  procedure  of  the  two 
individual  lateral  channels  as  developed  by  Arrow  [Ref.  6].  The 


design  was  performed  using  classical  techniques  and  involved  the 
uncoupled  yaw  and  roll  channels  for  the  elliptical  and  circular 
airframes  respectively.  The  resulted  continuous  open  loop 
designs  were  analyzed  in  terms  of  their  transient  and  frequency 
responses  and  found  to  be  in  accordance  with  the  desired 
requirements  specified  in  [Ref.  6].  Then,  applying  analog-to- 
digital  conversion,  the  corresponding  discrete  lateral  autopilots 
were  obtained  and  analyzed. 

The  second  part  utilized  modern  control  design  methods  to  the 
already  discussed  discrete  single-input  lateral  autopilots.  This 
allowed  comparison  with  the  preceding  classical  design,  and  more 
importantly  established  a  technique  to  extend  some  of  the  results 
to  the  more  general  multivariable  case.  At  first,  assuming 
availability  of  all  states  for  feedback  purposes,  application  of 
the  Ackermann  formula  led  to  a  discrete  state-feedback  designed 
autopilot.  Next,  since  the  state  vector  of  the  state- feedback 
model  is  not  usually  accessible  to  direct  measurement,  an  estimator 
was  introduced  as  an  additional  dynamic  design  in  order  to 
implement  control  to  the  original  system.  The  state-feedback  and 
estimator  designed  autopilots  were  analyzed  for  both  lateral 
channels  and  found  to  have  satisfactory  responses. 

Finally,  coupling  the  discrete  pitch  and  roll  channel 
autopilots,  the  state-feedback  and  estimator  designs  were 
obtained  and  proved  to  be  robust. 

The  analysis  in  all  the  above  cases  was  performed  using  the 
existed  at  Naval  Postgraduate  School  OPTSYS  and  ORACLS  Fo.  tran 


II .  CLASSICAL  DESIGN  AND  ANALYSIS  OF  LINEAR  UNCOUPLED 

LATERAL  AUTOPILOTS 


A.  GENERAL 

The  initial  phase  in  the  design  of  lateral  CBTT  autopilots 
involved  the  design  and  analysis  of  the  individual  uncoupled 
lateral  charnels  (i.e.,  yaw  and  roll)  with  prescribed 
relationships  between  speeds-of-response.  These  relationships 
when  coupled  with  the  corresponding  ones  of  the  longitudinal 
uncoupled  channel  (i.e.,  pitch)  would  meet  the  requirements  of 
the  overall  CBTT  autopilot. 

The  uncoupled  autopilot  design  method  was  classical  and  used 
a  combination  of  frequency  response  and  root  locus  techniques. 
Utilization  of  this  particular  design  method  led  to  the 
achievement  of  practical  bandwidths  (i.e.,  sufficient  high 
frequency  attenuation),  and  in  turn  provided  the  range  of  requir 
missile  body  angular  rates  and  control  motions.  In  addition,  th 
resulting  design  minimized  the  influence  of  aerodynamic 
variations  on  desired  responses.  The  application  of  the 
uncoupled  channels  to  the  whole  CBTT  autopilot  was  accomplished 
by  an  appropriate  choice  of  the  relative  time  constants  of  the 
individual  channels.  In  order  to  achieve  the  desired  maneuver 
plane  acceleration  the  roll  channel  was  designed  to  have  a  time 
constant  of  0.5  seconds.  The  yaw  uncoupled  channel,  which 
follows  the  roll  motion  to  produce  the  required  coordination 
(i.e.,  minimization  of  sideslip  angle),  is  designed  to  have  a 


more  rapid  response  with  a  time  constant  of  0.39  seconds  for  the 
circular  airframe.  The  detailed  requirements  for  the  classical 
design  of  the  uncoup’ ed  yaw  and  roll  channel  autopilots  are 
presented  in  Appendix  A. 

A  fixed  flight  condition  (i.e.,  constant  Mach  number  and 
altitude)  was  selected  for  this  preliminary  performance  study. 
Fixed  flight  conditions  are  typically  used  in  autopilot  designs 
to  identify  and  cure  critical  areas  of  concern.  When  autopilot 
requirements  are  satisfied  at  fixed  flight  conditions,  then  areas 
of  concern  caused  by  varying  them  are  addressed.  The  selected 
flight  condition  of  60000  feet  altitude  and  Mach  number  3.95 
provided  sufficient  dynamic  pressure,  so  that  the  missile 
maneuvers  resulted  in  large  enough  angles-of-attack  to  exercise 
sideslip  control.  Aerodynamic  data  for  this  particular  flight 
condition  are  provided  in  Appendix  A. 

The  aerodynamic  models  developed  for  stability  studies  in  the 
frequency  domain  were  linearized  about  a  trim  angle-of-attack  for 
both  lateral  channels.  The  following  three  assumptions  were 
made  : 


1.  The  plane  of  Figure  2.1  was  the  maneuver  plane. 

2.  The  missile  was  trimmed  in  pitch  (i.e.,  M  y  =  0,  at  fixed 
values  of  X  ,  q,  and  Cp. 

Rather  than  use  the  assumption  that  the  missile  roll  rate  is 
approximately  zero  as  it  is  normally  done  for  the  roll  stabilized 
3TT  control,  the  following  assumption  was  made  for  BTT: 


Axis 


K' 


¥ 
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3.  The  missile  roll  rate  was  constant. 

Linearized  aerodynamic  derivatives  are  given  in  Appendix  C. 

In  this  chapter  the  analysis  of  both  continuous  and  discrete 
uncoupled  lateral  channels  of  a  BTT  autopilot  was  based  on  the 
transient  and  frequency  responses  of  maneuver  plane  accelerations, 
body  angular  rates  and  tail  incidence  angles.  A  general  block 
diagram  of  a  BTT  autopilot  with  all  its  channels  is  shown  in 
Figure  2.2.  Inertial  acceleration  command  were  applied  in  polar 
coordinates  (i.e.,  magnitude  of  the  command  applied  to  pitch 
and  the  direction  4:^  to  roll  autopilot).  The  yaw  autopilot  was 
slaved  to  the  roll  autopilot  in  order  to  minimize  the  yaw  and 
roll  motions.  Achieved  maneuver  plane  accelerations  in 
rectangular  coordinates  (i.e.,  and  )  were  determined  by 

resolving  achieved  body-fixed  accelerations  (i.e.,  /i ^  and  /j,.  ) 
through  missile  roll  rate  (i.e.,  Euler  angles  0  and  were 
assumed  to  be  sufficiently  small). 

B.  AIRFRAME  CONFIGURATIONS 

The  two  airframe  configurations  studied  in  this  work  were 
taken  from  [Ref.  6]  and  are  shown  in  Figure  2.3  and  Figure  2.4. 
Although  the  configuration  in  Figure  2.3  reveals  a  body  of 
circular  cross  section  and  that  of  Figure  2.4  an  elliptical  one, 
both  airframes  have  the  same  cross  sectional  area  distribution. 

In  specific,  the  circular  cross  sectional  body  has  a  closure 
ratio  Abase/Am ax  of  O.69  with  Am ax  occurring  at  585  missile  body, 
whereas  the  elliptical  airframe  has  a  3:1  cross  section. 
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Figure  2,2  Bank 


gure  2.3  Model  of  Circul 


Both  airframe  configurations  are  tail-controlled  using  four 
identical  control  surfaces  which  are  located  flush  with  the  body 
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base  with  a  +30°  dihedral.  In  the  case  of  the  elliptical  body, 
the  hinge  line  was  skewed  such  that  a  10°  control  deflection 
measured  at  the  body-tail  juncture  had  a  resultant  7.04°  surface 
deflection.  Thus,  the  aerodynamic  control  effectiveness  in  terms 
of  deflection  measured  at  the  body-tail  is  lower  for  the 
elliptical  airframe  although  it  is  nearly  the  same  in  terms  of 
resultant  surface  deflection. 

The  total  span  of  the  mono-wings  is  the  same  for  each 
configuration,  which  results  in  larger  wing  area  for  the  circular 
airframe.  The  wing  area  and  span  for  the  circular  airframe  were 
chosen  as  typical  of  current  maneuvering  missiles.  The  wing  for 
the  elliptical  concept  was  determined  by  projecting  the 
elliptical  body  on  the  circular  body-wing  planform.  The  resultant 
when  exposed  wing  planform  became  the  wing  for  the  elliptical 
body . 

Comparison  of  the  elliptical  airframe  with  the  corresponding 
circular  indicates  the  following: 

About  30^  more  normal  force  that  is  nearly  independent  of 
angle-of-attack  can  be  achieved  at  supersonic  speeds. 

2.  Values  of  longitudinal  stability  oarameter  C,~  are  more 
positive,  and  with  more  pronounced  nonlinearities  in 
pitching  moment  at  subsonic  speeds. 

3.  Levels  of  directional  stability  are  increased  and  more 
compatible  with  levels  of  longitudinal  st'bility. 


31 


I 


4.  More  yaw  control  is  available  although  suitable  locations 
for  tails  on  the  body  are  more  limited  because  of  the 


geometry  of  the  elliptical  airframe. 

The  two  airframe  configurations  were  sized  to  provide 
realistic  geometric  and  mass  properties.  The  details  are 
presented  in  Appendix  D. 

C.  UNCOUPLED  YAW  CHANNEL  AUTOPILOT  FOR  ELLIPTICAL  AIRFRAME 

The  purpose  of  the  uncoupled  yaw  channel  autopilot  of  a  CBTT 
missile  is  to  minimize  the  sideslip  angle  (  j*  ) ,  or  provide 
coordinate  motion  between  the  yaw  and  roll  channels.  The  easiest 
way  to  accomplish  this  is  by  designing  the  uncoupled  yaw  channel 
(i.e.,  roll  and  pitch  dynamic  effects  neglected)  as  a  regulator 
(i.e.,  no  guidance  command  and  with  rate  and  acceleration 
feedback)  to  help  minimize  the  sideslip  angle. 

A  block  diagram  of  the  uncoupled  yaw  channel  is  shown  in 
Figure  2.5.  In  this  diagram  both  the  aerodynamic  model  and  yaw 
control  law  are  involved.  The  normal  acceleration  (  )  is  not 

used  to  command  the  CBTT  autopilot.  Instead,  it  is  used  for  the 
design  and  analysis  of  the  uncoupled  channel.  The  command  used  by 
the  coupled  system  is  shown  in  dashed  lines  and  is  a  yaw  angular 
rate  command  (r^).  The  yaw  control  law  shown  in  Figure  2.6  [Ref. 
6]  is  governed  by  missile  body  angular  rate  (r)  and  yaw  normal 
acceleration  (li.,,).  At  the  flight  condition  of  interest  (i.e., 

60  kft  altitude  and  Mach  number  3.95)  the  yaw  control  law 
determines  the  required  command  ( 
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)  to  an  actuator  which  is 


-  _  (S7.5) 


(II .C. 1-6 ) 


K  = 


(65+ 

V 


(II. C. 1-7) 


Substituting  the  values  of  aerodynamic  data 
linearized  aerodynamic  derivatives  (Table  VIII),  and 
and  mass  properties  (Table  IX),  equations  (II. C. 1-3) 
( II . C . 1 -7 )  become : 

A  =  -0.1351 
B  =  0.0436 
C  =  -17.2748 
E  =  34.5495 
K  =  0.4823 


(Table  VII), 

geometric 

through 


(II  .C.  1-8) 
(II. C.  1-9) 
(II. C. 1-10) 
(II. C. 1-1 1 ) 
(II .C. 1-12) 


Introducing  the  above  equations  (II. C. 1-8)  through 
(II. C. 1-12)  to  (II. C. 1-1)  and  (II. C. 1-2),  the  aerodynamic 
transfer  functions  of  the  uncoupled  yaw  channel  for  the  circular 
airframe  and  zero  angle-of-attack  can  be  obtained. 

a.  Transfer  Function  of  Yaw  Angular  Rate: 


I  ^  0.  0  US 


S  +i 


(II. C. 1-13) 


b.  Transfer  Function  of  Yaw  Normal  Acceleration: 

>  ^  I  f  i) 


(II .C. 1-14 ) 


Rearranging  the  above  transfer  functions  in  terms  of  the 
variable  pairs  (r,  -y  )  o^.)  respectively,  they  become: 


<rsUo.o4‘-/4ZS  -1S.^Si5r.--5c.^6'oyS  -^.7-4^Uly  (II. C.  1-15) 

lyS  ?.4^^SSy/-0.O01S  ryS-  IH.  6SO^  Cy  (II. c. 1-16) 


Applying  inverse  Laplace  transformation,  the  following 


set  of  linear  differential  equations  can  be  obtained: 


(II  .C. 1-17) 


Iv  f  L,C4‘^4rly-1!i..^S'  15  hy  T  ly  *  I'-CClSC-y  -  6SC1  c  y  Qj  ^  1-18) 


Both  equations  (II. C. 1-17)  and  (II.C.1-18)  form  a  second- 
order  system  of  linear  differential  equations  in  which  the 
forcing  function  involves  derivative  terms.  Using  rules  of 
state-space  representation  [Ref.  8]  the  following  equations  are 
obtained : 


X,  =  '  J’6.  c-y 

Xj  -  1^,-^Slb  0.04 ‘^4  Xj  “  U,  Cy 

r  c  -aii34  oy 

=  IS. 0^-74  Z^-^S'G.^TbbSc-y 

r  =  -  ^6,  66  Cy 


(II .C. 1-19) 
(II. C. 1-20) 
(II .C. 1-21 ) 
(II. C. 1-22) 
(II .C. 1-23) 
(II. C. 1-24) 


t-,  -  -^1  i  2 .  A62o  fly 


(II. C. 1-25) 


2. 


ions  of  Yaw  Control  Law  and  Actuato 


Equat 

a.  Acceleration  Compensator  Equation 

The  acceleration  compensator  equation  obtained  from 


Figure  2.6  is: 


Rearranging  and  applying  inverse  Laplace  transforma¬ 
tion,  (II. C. 2-1)  turns  into  the  following  linear  differential 
equation : 

y  =  -  5/  -t- 1 .  S-y  "o  i'l  y  -  1,  59/3  -)  (I  I .  C  .  2-2  ) 

Substituting  equation  (II. C. 1-25)  into  (II. C. 2-2), 
the  last  becomes: 

y  =  1, 59/5 Z  -  5 y  +  3.  y  554-/y  -  ;.5<9r5  ily^  (II .  C .  2-3  ) 

b.  Rate  Compensator  Equation 

The  rate  compensator  equation  also  obtained  from 


Figure  2.6  is: 

4,55  (’-Lin)  (II.c.2-!i) 

'T,  =  - - -  <■>'<■^1 

Rearranging  and  applying  inverse  Laplace  transforma¬ 
tion  to  (II. C. 2-5),  it  turns  into: 

=  L.  *L,4c5  r  f  4. S5y  f  4-, 55  ^  ( ll  .C.  2-5  ) 

Substituting  equations  (II. C. 1-23),  (II. C. 1-24)  and 
(iI.C.2-3)  into  the  above,  it  becomes: 

-  ,  -  -r  o5  X,  5x46,  (r+T’Z  -I  /4<5-/.15  SCM  :v'v..T/TrH-  (II. C. 2-6) 


o.  Actuator  Equation 

The  actuator  equation  obtained  from  Figure  2.5  is; 

1 

■  - ; -  (II. C. 2-7) 

Rearranging  and  applying  inverse  Laplace  transforma¬ 
tion,  the  last  equation  turns  into: 

SyrigB.  (II.C.2-8) 

3 .  Design  Approach  and  Analysis  of  Continuous  System 
Utilizing  state-space  representation,  the  equations 
(II. C. 1-19)  through  (II. C. 1-22),  (II. C. 2-2)  (II. C. 2-6)  and 
(II.C.2-8)  can  be  modeled  in  a  seventh-order  system  of  the  form 
x=Fx+Gu.  The  continuous  plant  system  and  input  matrixes  F  and  G 
in  Table  I,  and  the  state  vector  is: 


(II. C. 3-1  ) 

state  variables  are: 
yaw  angular  rates 
yaw  normal  accelerations 
output  of  acceleration  compensator 
input  command  in  the  actuator 


are  shown 

y 

'  C  V 

I 

Lv, 

i-  ' 

where  the 
,  X2  : 

z  1  ,  Z2  : 
y 

W  V  • 

ic 


yaw  tail  incidence 


PLANT  SYSTLM  AND  INPUT  MATRICliS;  UNCOUPLLD  YAW 
CHANNEL  AUTOPILOT;  CLASSICAL  DESIGN;  CONTINUOUS  OPEN 


Executing  the  OPTSYS  program,  using  an  input  step 
function  which  represents  "1  gee  command"  at  zero  trim  angle-of- 
attack,  the  pole-zero  and  time  and  frequency  response  plots  are 
obtained . 

The  pole-zero  plot  of  Figure  2.7  indicates  that  the 


continuous  open  loop  system  is  stable,  since  the  s-plane  poles 
are : 


S,-  -  |74.S7<S 

(II. C. 3-2) 

5^  =  '6.  !04^  +-J 

(II. C. 3-3) 

lC.^5v6 

(II. C. 3-4) 

S4  ■=  -C,C2i74-7S  4-,CSvI‘v 

(II. C. 3-5) 

Si  r-t;.C2l74-75-j4CS9l^ 

(II. C. 3-6) 

(II. C. 3-7) 

S;  :  9’c  -j  ^  4 

(II. C. 3-8) 

The  time  response  plots  of  the  yaw  normal  acceleration,  angular 
rate  and  tail  incidence  are  shown  in  Figures  2.8  through  2.10. 

In  particular  the  yaw  normal  acceleration  time  response  plot  has 
a  0.39  seconds  time  constant,  1%  overshoot  and  a  steady-state 
error  of  0.018.  These  results  are  in  accordance  with  the 
requirements  referred  in  Appendix  A,  that  is  a  time  constant  of 
0.4  seconds,  overshoot  less  than  10X  and  steady-state  error  not 
necessarily  equal  to  zero.  AIL  the  above  three  time  response 
plots  are  identical  with  those  presented  in  [Ref.  6]. 

Figures  2.11  through  2.16  show  the  frequency  response 
plots  of  the  yaw  normal  acceleration,  angular  rate  and  tail 
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incidence,  from  which  the  phase  crossover  frequencies  and  gain 
margins  of  Table  II  can  be  obtained.  The  positive  gain  and  phase 
margins  of  the  open  loop  system  ensure  the  relative  stability  of 
the  closed  loop  (controlled)  system. 

4 .  Design  Approach  and  Analysis  of  Discrete  System 

Utilizing  analog- to-digital  conversion  by  the  aid  of 
ORACLS  program  and  for  a  sample  period  of  0.0125  seconds,  a 
seventh-order  discrete  system  of  the  form  x(k+1  )  =  Ax(k)+Bu(k)  is 
obtained.  The  discrete  pant  system  and  input  matrices  A  and  B 
are  shown  in  Table  III. 

The  pole-zero  plot  of  Figure  2,1?  indicates  that  the 
discrete  open  loop  system  is  also  stable,  since  the  z-plane 
poles  are: 


z,  =  C.\\^C7€ 

1 

(II. C. 4-1 ) 

Z2.  -  9l6554-tJ 

(II. C. 4-2) 

Z5  --  C. ‘-7 1  & ^ S4 - j a  \ZXB6l 

(II. C. 4-3) 

^.99^411  f  J  O.Oii6V6a 

(II. C. 4-4) 

Z3-  --  0  9984-11  -j 

(II. C. 4-5) 

--  C.9t55CSi-jO.O^C  I^48 

(II. C. 4-6) 

Zj  --  L.  94  S56  S -J  C  C^cT  §46 

(II.C.4-7) 

The  time  response  plots  of  the  yaw  normal  acceleration, 
angular  rate  and  tail  incidence  for  the  discrete  uncoupled  yaw 
channel  are  presented  in  Figures  2.18  through  2.20.  A  close 
observation  of  these  plots  indicates  that  they  are  identical  with 
those  of  the  continuous  classical  system  found  in  the  previous 
section . 
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Figure  2.7 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel 
Classical  Design  Continuous  Open  Loop 
Elliptical  Airframe 


Autopilot ; 
System ; 
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Figure  2. 


9  Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Elliptical  Airframe 


JCED  AI  GOV  Lf^rj'.U  tJ  I  fxi’crjst 


MAGNITUDE  -  DECIBELS 


FREQUENCY  -  RADIANS  PER  SECOND 


Yaw  Angular  Rate -Gain  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Air f  rarae 
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FREQUENCY  -  RADIANS  PER  SECOND 


Figure  2.14  Yaw  Angular  Rate-Phase  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2,15  Yaw  Tail  Incidence-Gain  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airframe 
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FREQUENCY  -  RADIANS  PER  SECOND 


Yaw  Tail  Incidence-Phase  vs  Frequency;  Uncoupled 
Yaw  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Elliptical 
Airf  name 


PKASli  CROSSOVER  I'RliQUENCIRS  AND  GAIN  MARGINS; 
OllPl.P.n  YAW  CilANNRL  AUTOPILOT;  CLASSICAL  DLSIGN; 


D.  UNCOUPLED  POLL  CHANNEL  AUTOPILOT  FOR  CIRCULAR  AIRFRAME 

The  uncoupled  roll  channel  autopilot  of  a  CBTT  missile  is 
commanded  to  roll  the  missile  so  as  to  put  the  preferred  maneuver 
direction  in  the  direction  of  the  guidance  acceleration  command. 
The  desired  maneuver  plane  acceleration  should  be  attained  as 
rapidly  as  the  achieved  body-fixed  pitch  acceleration.  To 
accomplish  this,  the  uncoupled  roll  channel  autopilot  (i.e.,  yaw 
and  roll  dynamic  effects  neglected)  was  designed  to  have  the  roll 
angle  time  constant  equal  to  the  time  constant  of  the  normal 
acceleration  achieved  by  the  uncoupled  pitch  channel  autopilot. 

A  block  diagram  of  the  uncoupled  roll  channel  is  shown  in 
Figure  2.21.  In  this  diagram  both  the  aerodynamic  model  and  roll 
control  law  are  involved.  The  roll  control  law  shown  in  Figure 
2.22  is  commanded  by  roll  angle  (|)q)  and  governed  by  roll  angular 
rate  (p)  and  roll  angle  (5j!). 


The  design  and  analysis  of  the  uncoupled  roll  channel 
autopilot  was  performed  in  this  section  for  the  stable  at  zero 
angle-of-attack  circular  airframe. 

1 .  Transfer  Functions  of  Aerodynamic  Model 

The  aerodynamic  transfer  functions  of  the  uncoupled  roll 
channel  autopilot  obtained  from  Figure  2.21  are;- 

a.  Transfer  Function  of  Roll  Angular  Rate: 

_  i  I  O  ll  r.- 

_  - - — >_ 


(II. D. 1-1 ) 


PLANT  SYSTLM  AND  INPUT  MATRICES;  UNCOUPLED  YAW 
CHANNEL  AUTOPILOT;  CLASSICAL  DESIGN;  DISCRETE  OPEN 
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Figure  2.17  Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot 
Classical  Design;  Discrete  Open  Loop  System; 
Elliptical  Airframe 
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Yaw  Angular  Rate  vs  Time;  Uncoupled  Y 
Channel  Autopilot;  Classical  Design; 
Discrete  Open  Loop  System;  Elliptical 
Airframe 
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Figure  2.20  Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Classical  Design; 
Discrete  Open  Loop  System;  Elliptical 
Airframe 


b. 


Transfer  Function  of  Roll  Angle: 


_ 

^  = 


1 


(II  .D.  1-2  ) 


Substituting  the  values  of  aerodynamic  data  (Table  VII), 
linearized  aerodynamic  derivatives  (Table  VIII)  and  geometric  and 
mass  properties  (Table  IX),  equation  (II. D.  1-1)  becomes 

ps  =  8.03462  (II. D.  1-3) 

Applying  inverse  Laplace  transformation  to  equations 
(II. D. 1-2)  and  (II. D. 1-3),  the  following  set  of  linear 
differential  equations  is  obtained: 

p  =  8.03462  (II  .D,  1-4) 

9  =  p  (II. D. 1-5) 

2 .  Equations  of  Roll  Control  Law  and  Actuator 
a.  Roll  Angle  Compensator  Equation 

The  roll  angle  compensator  equation  obtained  from 

Figure  2.22  is: 

-) 

4 


X 


(II. D. 2-1  ) 


Rearranging  and  applying  inverse  Laplace 
transformation  (II. D. 2-1)  turns  into  the  following  linear 
differential  equation: 

X  =  -17. 6^?-  8X  +  17.6C{:^  (II. D. 2-2) 

b.  Rate  Compensator  Equation 

The  rate  compensator  equation  also  obtained  from 


Figure  2.22  is: 

Y_  'iisccb  (  — 


CX'P) 


(II. D. 2-3) 


I  t  ^  i ;  .  f  j  ^  1  1  'v  ^ ‘  ^  IV  ‘  I  j  ) » I  Cl  o  d  J  -S 


Fij;ure  2.21  Uncoupled  Roll  Channel 


Rearranging  and  applying  inverse  Laplace  transforma¬ 


tion  to  (II. D. 2-3),  it  turns  into: 

Y  +  5Y  =  0.0502X  -  0.0502p  +  0.75325X  -  0.75325p  (II. D. 2-4) 

Substituting  equations  (II. D. 1-4)  and  (II. D. 2-2)  into 
(II. D. 2-4),  the  last  becomes: 

Y=-0.7532p-0. 88352CH0. 35 165X-5Y-0. 40334  Og+0.88352'?'(,  (II.D.2-5) 

G.  Pseudo-Differential  Equation 

The  Dseudo-dif ferential  eauation  obtainpd  from  Fivura 


2.22  is: 

:.  Cl  50??  5 


P 


(II. D. 2-6) 


Rearranging  and  applying  inverse  Laplace 
transformation  to  (II. D. 2-6),  it  turns  into: 

X^  =  -6Xi  +  0.078l98p  (II. D. 2-7) 

Substituting  equation  (II. D. 1-4)  into  (II. D. 2-7),  the 
last  becomes: 


X^  =  -  6Xt  +  0.628291 


(II .D. 2-8) 


d.  Equation  of  Actuator  Compensator 

The  equation  of  actuator  compensator  obtained  from 
Figure  2.22  is: 

S 


lie 


fl 


s 


IS 


f  i 


c  r-Xi) 


(II .D.2-9  ) 


Rearranging  and  applying  inverse  Laplace 
tion  to  (II. D.2-9),  it  turns  into; 
o.  =-15^'^  +0.  1  3636Y-0.  1  3636Xi  +  15Y-15Xi 


transforma- 

(II . D. 2-10) 


Substituting  equations  (II. D. 2-5)  and  (II. D, 2-8)  into 


(II. D. 2-10),  the  last  becomes: 

-  -  C.  IC^JI  p  -C  UC4Scf  i-C  C414SX<r  14:  Z\SZ  r-l«4;  '  1 5  O.  mc  J-Cg  f  CM  M 

(III  .D. 2-11 ) 


e.  Actuator  Equation 

The  actuator  equation  obtained  from  Figure  2.21  is: 


^r/.5 


_i _ 

lffc.4 


5  4-1 


(II. D. 2-12) 


Rearranging  and  applying  inverse  Laplace 
transformation,  the  last  equation  turns  into: 


=  10795.32&j^  -188.4 


(II. D. 2-13) 


3 •  Design  Approach  and  Ana'ysis  of  Continuous  System 
Utilizing  state-space  representation,  the  equations 


(II. D. 1-4),  (II. D. 1-5),  (II. D. 2-2),  (II. D. 2-5),  (II. D. 2-8), 

(II. D. 2-11),  and  (II. D. 2-13)  can  be  modeled  in  a  seventh-order 
system  of  the  form  x=Fx+Gu.  The  continuous  plant  system  and 
input  matrices  F  and  G  are  shown  in  Table  IV,  and  the  state 


vector  is: 


XcK)  = 


? 

X 

Y 

X. 


(II. D. 3-1 ) 


where  the  state  variables  are: 


roll  angular  rate 


(j)  :  roll  angle 

X  :  output  of  roll  angle  compensator 

Y  :  output  or  rate  compensator 

X  :  output  of  pseudo-differentiator 

0|^  :  input  command  in  the  actuator 

:  roll  tail  incidence 

Executing  the  OPTSYS  program,  using  an  input  step 
function  which  represents  "1  gee  command”  at  zero  trim  angle-of- 
attack,  the  pole-zero  and  time  and  frequency  response  plots  are 
obtained . 

The  pole-zero  plot  of  Figure  2.23  indicates  that  the 
continuous  open  loop  system  is  stable,  since  the  s-plane  poles 


f  =  -174.785  (roll  angular  rate) 
5,  =  -9.25097+028.4098  (roll  angle) 


(II. D. 3-2) 
(II. D. 3-3) 


^5  =  -9 . 25097- j28 . 4098  (output  of  roll  angle  compensator) 

(II. D. 3-4) 

=  -2 . 46608+ j2 . 7 1 1 52  (output  of  rate  compensator  network) 

(II. D. 3-5) 

=  -2 . 46608- j2 . 7 1 1 52  (output  of  pseudo-differentiator) 

(II. D. 3-6) 


=  -8.98209  (input  command  in  the  actuator) 
i  j.  =  -5.19932  (roll  tail  incidence) 


(II. D. 3-7) 
(II. D. 3-8) 


The  time  response  plots  of  the  roll  angle,  angular  rate 
and  tail  incidence  are  shown  in  Figure  2.24  through  2.26.  In 
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particular,  the  roll  angle  time  response  plot  has  a  0.55  seconds 
time  constant,  3^  overshoot  and  a  steady-state  equal  to  zero. 
These  results  are  in  accordance  with  the  requirements  referred  in 
Appendix  A,  that  is  a  time  constant  of  0.5  seconds,  overshoot 
less  than  10^  and  zero  steady-state  roll  angle  error.  All  the 
above  three  time  response  plots  are  identical  with  those 
presented  in  [Ref.  6]. 

Figures  2.27  through  2.32  show  the  frequency  response 
plots  of  the  roll  angle,  angular  rate  and  tail  incidence,  from 
which  the  phase  crossover  frequencies  and  gain  margins  of  Table 
V  can  be  obtained.  The  positive  gain  and  phase  margins  of  the 
open  loop  system  ensure  the  relative  stability  of  the  closed  loop 


(controlled)  system. 
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Figure  2.26  Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Continuous  Open  Loop 
System;  Circular  Airframe 
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Figure  2.27  Roll  Angle-Gain  vs  Frequency;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Continuous  Open  Loop 
System;  Circular  Airframe 


FREQUENCY  -  RADIANS  PER  SECOND 


Figure  2.28  Roll  Angle-Phase  vs  Frequency;  Uncoupled  Roll 

Channel  Autopilot;  Classical  Design;  Continuous 
Open  Loop  System;  Circular  Airframe 


FREqUENCY  -  RADIANS  PER  SECOND 


Roll  Angular  Rate -Gain  vs  Frequency;  Uncoupled  Roll 
Channel  Autopilot;  Classical  Design;  Continuous  Open 
Loop  System  Circular  Airframe 


FREQUENCY  -  RADIANS  PER  SECOND 


Roll  Angular  Rate-Phase  vs  Frequency;  Uncoupled  Roll 
Channel  Autopilot;  Classical  Design;  Continuous  Open 
Loop  System;  Circular  Airframe 


MAGNITUDE  -  DECTDELS 


FREQUENCY  -  RADIANS  PER  SECOND 


Roll  Tail  Incidence-Phase  vs  Frequency;  Uncoupled 
Roll  Channel  Autopilot;  Classical  Design; 
Continuous  Open  Loop  System;  Circular  Airframe 


TABLE  V 


PHASE  CROSSOVER  FREQUENCIES  AND  GAIN  MARGINS;  UNCOUPLED  ROLL 
CHANNEL  AUTOPILOT;  CLASSICAL  DESIGN  CONTINUOUS  OPEN  LOOP  SYSTEM; 

CIRCULAR  AIRFRAME 


PHASE  CROSSOVER 

GAIN 

FREQUENCY  (rad/sec) 

MARGIN  ( 

ROLL 

ANGLE  (cp) 

10.0461 

60.2165 

ROLL 

ANGULAR  RATE  (p) 

10.0461 

60.2165 

ROLL  TAIL  INCIDENCE  ) 


4 .  Design  Approach  and  Analysis  of  Discrete  System 

Utilizing  analog- to-digital  conversion  by  the  aid  of 
ORACLS  program  and  for  a  sample  period  of  0.0125  seconds,  a 
seventh-order  discrete  system  of  the  form  x(K+1 )=Ax(k)+  Bu(k)  is 
obtained.  The  discrete  plant  system  and  input  matrices  A  and  B 
are  shown  in  Table  VI. 

The  pole-zero  plot  of  Figure  2.33  indicates  that  the 
discrete  open  loop  system  is  also  stable,  since  the  z-plane  poles 
are : 

-  0.0992805  (roll  angular  rate)  (II. D. 4-1) 

-  0.8352l6+j0. 309736  (roll  angle)  (II. D. 4-2) 

2’?  =  0 . 8352 1  6- jO  .  309736  (output  of  roll  angle  compensator) 

(II. D. 4-3) 

24  =  0 . 969087+JO . 0328588  (output  of  rate  compensator  network) 

(II .D. 4-4  ) 

Z5  =  0.969087-10.0328588  (output  of  pseudo-differentiator) 

(II .D. 4-5 ) 

2^  =  0.893797  (input  command'  in  the  actuator)  (II. D. 4-6) 

-  0.937075  (roll  tail  incidence)  (II. D. 4-7) 

The  time  response  plots  of  the  roll  angle,  angular  rate 
and  tail  incidence  for  the  discrete  uncoupled  roll  channel  are 
presented  in  Figures  2.34  through  2.36.  A  close  observation  of 
these  plots  indicates  that  they  are  identical  with  those  of  the 
continuous  classical  system  found  in  the  previous  section. 
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Figure  2.35  Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Classical  Design;  Discrete  Open  Loop 
System;  Circular  Airframe 
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III.  MODERN  CONTROL  DESIGN  AND  ANALYSIS  OF  LINEAR 
UNCOUPLED  LATERAL  AUTOPILOTS 


A.  GENERAL 

The  task  of  this  chapter  is  the  design  and  analysis  of  the 
same  discrete  uncoupled  lateral  channel  autopilots  discussed  in 
the  previous  chapter,  using  different  techniques  which  are  based 
on  modern  control  formulation.  The  difference  in  the  two 
approaches  is  entirely  in  the  design  method  since  the  end  result, 
a  set  of  difference  equations  providing  control,  is  identical. 

Modern  control  theory  is  contrasted  with  the  classical 
control  theory  in  that  the  former  is  applicable  to  multi-input- 
multi-output  systems,  which  may  be  linear  or  nonlinear  time- 
invariant  or  time-varying,  while  the  latter  is  applicable  only  to 
linear  time-invariant  single-input-single-output  systems.  Also, 
modern  control  theory  is  essentially  a  time-domain  approach, 
while  the  conventional  classical  control  theory  is  a  complex 
frequency-domain  approach. 

System  design  in  classical  control  theory  is  based  on  trial- 
and-error  procedures  which,  in  general,  will  not  yield  optimal 
control  systems.  System  design  in  modern  control  theory,  on  the 
other  hand,  enables  the  design  of  optimal  control  systems  of 
great  complexity  and  good  accuracy  with  respect  to  given 
performance  indexes.  In  addition,  design  in  modern  control 
theory  can  be  carried  out  for  a  class  of  inputs  instead  of 
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specific  input  function,  such  as  the  impulse,  step  or  sinusoidal 
functions  and  can  also  include  initial  conditions. 

One  of  the  most  attractive  features  of  modern  control  design 
method  is  that  the  procedure  consists  of  two  independent  steps. 

One  step  assumes  that  all  the  system  states  are  available  for 
feedback  purposes.  In  general,  even  if  this  is  not  a  practical 
enough  assumption  since  it  needs  a  large  number  of  sensors,  it  is 
usually  adopted  in  order  to  accomplish  the  first  design  step, 
namely  the  control-law.  The  remaining  step  is  the  design  of  an 
estimator  which  estimates  the  entire  state  vector,  given 
measurements  of  portion  of  the  state  provided  by  the  system  output 
equation.  The  final  control  algorithm  consists  of  the  control- 
law  and  estimator  combined,  where  the  control-law  calculations 
are  based  ^n  the  estimated  states  rather  than  the  actual  states. 
This  substitution  is  reasonable  and  the  combined  design  can  give 
closed  loop  characteristics  which  are  unchanged  from  those 
assumed  in  designing  the  control-law  and  estimator  separately. 


B.  DISCRETE  STATE-FEEDBACK  DESIGN 

Considering  the  following  discrete  control  system: 

X(K  +  1)  =  Ax(k)  +  Bu(k)  (III. B. 1-1) 

Y(k)  =  Hx(k)  (III.B.  1-2) 

the  control-law  design  is  also  referred  to  as  state-feedback 
design  since  it  is  simply  the  feedback  of  a  linear  combination  of 
all  the  system  states,  that  is: 

u(k)  =  Fx(k)  (III.B. 1-3) 

where  F:  control-law  gain  vector 
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Thus,  the  characteristic  equation  of  the  controlled  (closed  loop) 
system  is: 


det  (zI-A+BF)  =  0  (III. B. 1-4) 

The  discrete  state-feedback  design,  providing  that  the  system  is 
controllable,  consists  then  of  finding  the  control-law  gain 
vector  F  so  that  the  roots  of  (III. B. 1-4)  are  in  desirable 
locations . 

A  program  logic  for  computing  the  control-law  gain  vector  via 
the  Ackermann's  formula  is  given  in  Appendix  E  [Ref.  9]. 

Utilizing  this  control  algorithm  a  Fortran  program  was  written 
(Appendix  F)  which  has  as  inputs  the  sample  period,  the  discrete 
plant  system  and  control  input  matrices,  the  s-plane  poles  and 
provides  as  output  the  control-law  gain  vector. 

1 •  Uncoupled  Yaw  Channel  for  Elliptical  Airframe 
a.  Control-Law  Gain  Vector 

Executing  the  Ackermann  Fortran  program  described  in 
Appendix  F  with  inputs: 

(1)  Sample  period  of  0.0125  seconds 

(2)  Discrete  plant  system  and  control  input  matrices  A  and  B 
of  Table  III 

(3)  S-plane  poles  defined  in  equations  (II. C. 3-2)  through 
(II. C. 3-8) 

the  following  control-law  gain  vector  for  the  elliptical  airframe 

of  the  uncoupled  yaw  channel  is  obtained: 

F=[-1.0195  -0.109  -0.3492  0.0393  -0.6152  32.4775  -31.0868] 

(III.B.  1-5) 


b.  Design  Approach  and  Analysis 


The  discrete  state-feedback  designed  yaw  autopilot 
can  be  found  by  introducing  the  control-law  gain  vector  of 
(II. B. 1-5)  into  the  original  system. 

The  pole-zero  plot  of  Figure  3.1  indicates  that  the 
discrete  closed  lop  system  is  stable,  since  the  z-plane  closed 
loop  poles  are: 

=  0.46461  (yaw  angular  rate)  (III. B. 1-6) 

^2,  =  0.898888  +  jO.  137039  (III. B. 1-7) 

■2’j-  =  0 . 898888- jO  .  1  37039  (yaw  normal  acceleration)  (III. B.  1-8) 

=  0.998423+j0. 051078  (III. B. 1-9) 

^5  =  0.998423- jO. 05 1 0789  (output  of  acceleration  compensator 

network)  (III. B. 1-10) 

-  0 . 962974+ jO , 035871 2  (input  command  in  the  actuator) 

(III.B.  1-11  ) 

=  0.962974-j0. 0358712  (yaw  tail  incidence)  (III.B. 1-12) 

The  time  response  plots  of  the  yaw  normal 
acceleration,  angular  rate  and  tail  incidence  are  presented  in 
Figures  3.2  through  3.4.  A  close  observation  of  the  above  pole- 
zero  and  time  response  plots  for  the  discrete  yaw  state-feedback 
design  indicates  that  they  are  identical  with  those  of  the 
discrete  classical  design  found  in  the  previous  chapter. 

G.  Simplified  Design 

The  discrete  state-feedback  designed  autopilot  of  the 
previous  section  can  be  simplified  by  reducing  the  returning  gain 
loops.  This  can  be  accomplished  by  placing  zeros  into  appropriate 
elements  of  the  control-law  gain  vector: 


F=[-1.0195  -0.109  -0.3492  0  -0.6152  32.4775  -31.0868] 

(III. B. 1-13) 

The  pole-zero  and  time  response  plots  of  the 
resulting  simplified  state-feedback  yaw  autopilot,  shown  in 
Figures  3.5  through  3-8,  do  not  present  significant  differences 
from  the  corresponding  plots  of  the  discrete  classical  design 
apart  from  the  overshoot  which  was  slightly  increased.  The 
discrete  closed  loop  system  is  again  stable,  since  the  z-plane 


closed  loop  poles  are: 

Z,  =  0.467984  1-14) 

=  0.891328+j0. 126019  (III. B. 1-15) 

Z5  =  0.891328-j0. 126019  (III.B.I-I6) 

<4.  =  0.998422+j0. 0510788  (III. B. 1-17) 

z-  =  0.998422-j0. 0510788  (III. B. 1.18) 

Z^  =  0.969075+00.0429695  (III. B. 1-19) 

=  0.969075-j0. 0429695  (III. B. 1-20) 

2 .  Uncoupled  Roll  Channel  for  Circular  Airframe 
a.  Control-Law  Gain  Vector 

Following  the  same  procedure  as  in  the  yaw  channel 
case  apart  from  the  use  of  discrete  matrices  A  and  B  from  Table 
III  and  s-plane  poles  from  (II. D. 3-2)  through  (II. D. 3-8),  the 
control-law  gain  vector  for  the  circular  airframe  roll  channel 


was  found  to  be: 

F=[0  0.0004  0  0.0178  -0.0179  -1.499  0.0023] 


(III .B. 2-1 ) 
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Figure  3-1  Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 

State-Feedback  Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 
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Figure  3. 


2  Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 

Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure  3. 


3  Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure  3,4  Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw  Channel 
Autopilot;  State -Feed back  Design;  Discrete 
Closed  Loop  System;  Elliptical  Airframe 
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Figure  3.5 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 
Simplified  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Elliptical  Airframe 
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Figure  3-6  Yaw  Normal  Acceleration  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  State-Feedback 
Design;  Discrete  Closed  Loop  System;  Elliptical 
Air  f  rame 
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ANGULAR  RATE.R 


b.  Design  Approach  and  Analysis 

The  discrete  state-feedback  designed  roll  autopilot 
can  be  found  by  introducing  the  control-law  gain  vector  of 
(III. B. 2-1)  into  the  original  system. 

The  pole-zero  plot  of  Figure  3.9  indicates  that  the 
discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 
loop  poles  are: 

-2,  =  0.0989714  (roll  angular  rate)  (III. B. 2-2) 

^2^  =  0.832909  +  j0. 315088  (roll  angle)  (III. B. 2-3) 

-  0 . 832909- jO . 31 5088  (output  of  roll  angle  compensator) 

(III. B. 2-4) 

=  0 . 969236+jO . 0333773  (output  of  rate  compensator  network) 

(III. B. 2-5) 

-  0 . 969236- jO . 0333773  (output  of  pseudo-differentiator) 

(III. B. 2-6) 

=  0.896583  (input  command  in  the  actuator)  (III. B. 2-7) 

=  0.938233  (roll  tail  incidence)  (III. B. 2-8) 

The  time  response  plots  of  the  roll  angle,  angular 
rate  and  tail  incidence  are  presented  in  Figures  3-10  through 
3.12.  A  close  observation  of  the  above  pole-zero  and  time 
response  plots  for  the  discrete  roll  state-feedback  design 
indicates  that  they  are  identical  with  those  of  the  discrete 
classical  design  found  in  the  previous  chapter. 

c.  Simplified  Design 

The  state-feedback  roll  designed  autopilot  can  be 


simplified  by  reducing  the  returning  gain  loops  as  follows: 

F=ro  0  0  0.178  -0.0179  -1.499  0]  (III. B  2-9) 
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Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot 
State-Feedback  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 
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Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure  3.12  Roll  Tail  Incidence  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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igure  3.14  Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot 
Simplified  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Roll  Angular  Rate  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Simplified  State-Feedback  Design; 
Discrete  Closed  Loop  System;  Circular  Airframe 
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Considering  the  same  as  in  the  state-feedback  case  discrete 
control  system,  a  prediction  estimator  defined  by  the  following 
equation  is  introduced: 

x(K+1)  =  Ax(k)  +  Bu(k)  +  K[y(k)  -  Hx(k)]  (III. C. 1-1) 

A 

where  x:  estimate  state  vector 
K:  estimate  gain  vector 

In  this  closed  loop  estimator,  shown  in  Figure  3.17,  the 
difference  between  the  measured  and  estimated  output  is  fed  back 
and  the  model  is  constantly  corrected  with  this  error  signal 

A 

which  is  defined  as  x=x-x.  The  difference  equation  describing 
the  behavior  of  the  error  is  obtained  by  subtracting  equation 
from  the  actual  plant  output  equation  (II. B. 1-2): 

x'(K  +  1)  =  [A  -  KH]  x(k)  (III. C. 1-2) 

Thus  the  characteristic  equation  of  the  controlled  (closed  loop) 
system  is: 

det  (zl  -  A  +  KH)  =  0  (III.C. 1-3) 

The  discrete  estimator  design,  providing  that  the  system  is 
observable,  consists  then  of  finding  the  estimator  gain  vector  K 
so  that  the  roots  of  (III.C. 1-3)  are  at  desirable  locations. 

The  estimator  gain  vector  can  be  obtained  again,  as  in  the 
case  of  the  state-feedback  design,  by  application  of  the 
Ackermann  program  of  Appendix  F  with  inputs  the  sample  period, 
the  transposes  of  the  discrete  plant  and  output  matrices,  and 
faster  s-plane  poles  from  those  of  the  continuous  system. 
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Executing  the  Ackermann  Fortran  program  (Appendix  F) 


with  inputs: 

(1)  Sample  period  of  0.0125  seconds 

(2)  The  transposes  of  the  discrete  yaw  system  plant  and  output 

matrices,  that  is  a"^  and  where: 


1 

0 

0 

0 

0 

0 

0 


(III .C. 1-4  ) 


(3)  S-plane  poles  slightly  faster  than  those  of  the  continuous 


open  loop  system,  that  is 
5,  =  -174.386 

5^  =  -6, 1 l42+j10.6396 
=  -6.  n42-j10.6396 
i.,.  =  -0.0227475  +  j4 .08919 
S5  =  -0.0227475-J4. 08919 
=  -2.9396+02.99929 
V  =  -2.9396-02.99929 


with  more  negative  real  parts: 

(III.C.  1-5) 
(III. C. 1-6) 
(III.C. 1-7) 
(III.C. 1-8) 
(III.C. 1-9) 
(III.C. 1-10) 
(III.C. 1-11) 


the  transpose  of  the  estimator  gain  vector  for  the  elliptical 
airframe  of  the  uncoupled  yaw  channel  is  calculated  as  output, 
from  which  the  following  K  is  obtained: 


K 


(III.C. 1-12) 


-0.0034 
0.0008 
-0.0003 
-0.0177 

0.0005 

-0.0017 

-0.0240 

b.  Design  Approach  and  Analysis 

The  discrete  estimator  designed  yaw  autopilot  can  be 
found  by  introducing  the  estimator  gain  vector  of  (III.C. 1-12) 
into  the  original  system. 

The  pole-zero  plot  of  Figure  3.18  indicates  that  the 
discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 
loop  poles  were  found  to  be: 

Z|  s  0.117115  (yaw  angular  rate)  (III.C. 1-13) 

^2.  -  0.918174+jO.  122926  (III.C. 1-14) 

■^5  =  0.9l8l74-j0. 122926  (yaw  normal  acceleration)  (III.C. 1-15) 

=  0.998399+JO. 0510958  (III.C. 1-16) 

25  =  0.998399- jO. 05 1 0958  (output  of  acceleration  compensator 

network)  (III.C. 1-17) 

=  0 . 963249+ jO . 0361 306  (input  command  in  the  actuator) 

(III.C.  1-18 ) 

-  0.963249-20.0361306  (yaw  tail  incidence)  (III.C. 1-19) 

The  time  response  plots  of  the  yaw  normal 
acceleration,  angular  rate  and  tail  incidence  are  presented  in 
Figures  3-19  through  3.21.  A  close  observation  of  the  above  pole 
zero  and  time  response  plots  of  the  discrete  yaw  estimator  design 
indicates  'hat  they  are  very  close  to  those  of  the  discrete 
classical  design  found  in  the  previous  chapter. 


Figure  3.13  Pole-Zero  Plot;  Uncoupled  Yaw  Channel  Autopilot; 

Estimator  Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 


Simplified  Design 

The  discrete  estimator  designed  autopilot  of  the 


previous  section  can  be  simplified  by  reducing  the  returning  gain 
loops.  This  can  be  accomplished  by  placing  zeros  into 
appropriate  elements  of  the  estimator  gain  vector: 


K 


-0.0034 

0 

0 

-0.0177 

0 

-0.0017 

-0.0240  j 


(III.C. 1-20) 


The  pole-zero  and  time  response  plots  of  the 
resulting  simplified  estimator  yaw  autopilot,  shown  in  Figures 
3.22  through  3.25,  do  not  present  significant  differences  from 
the  corresponding  plots  of  the  discrete  classical  design,  apart 
from  the  overshoot  which  was  slightly  increased.  The  discrete 
closed  loop  system  is  again  stable,  since  the  z-plane  closed  loop 
poles  are: 


=0.117117 

=  0. 9l8434+jO.  122898 
Zs  =  0. 9l8434-jO. 122893 
-^4  =  0.998422  +  j0. 0510954 
Z5  =  0.998422-j0. 0510954 
=  0.962965  +  j0. 0357874 
Zj  =  0.962965-j0. 0357874 


(III.C. 1-21 ) 
(III.C. 1-22) 
(III.C. 1-23) 
(III.C. 1-24) 
(III.C. 1-25) 
(III.C. 1-26) 
(III.C. 1-27) 


2 .  Uncoupled  Roll  Channel  for  Circular  Airframe 
a.  Estimator  Gain  Vector 

Following  the  same  procedure  as  in  the  yaw  channel 
case  apart  from  the  use  of: 

( 1 )  The  transpose  of  the  discrete  roll  system  plant  matrix 

(2)  S-Plane  poles  slightly  faster  than  those  of  the  continuous 


open  loop  roll  system,  that  is: 


= 

-184.795 

(III. C. 2-1 ) 

= 

-9.26097+j28.4098 

(III. C. 2-2) 

= 

-9.26097-j28.4098 

(III. C. 2-3) 

.54. 

= 

-2.47608+j2.71 152 

(III. C. 2-4) 

^■5 

= 

-2.47608-j2.71 152 

(III. C. 2-5) 

= 

-8.99209 

(III. C. 2-6) 

= 

-5.20032 

(III. C. 2-7) 

the 

estimator  gain  vector  for  the  circular  airframe 

oil  channel 

is  obtained: 


K 


0.004 

0 

0 

-0.0002 

0.0003 

-0.0492  J 

b.  Design  Approach  and  Analysis 

The  discrete  estimator  designed 


(III. C. 2-8) 


roll  autopilot  can  be 


found  by  introducing  the  estimator  gain  vector  of  (III. C. 2-3) 
into  the  original  system. 
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Figure  3.20  Yaw  Angular  Rate  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 


Pole-Zero  Plot;  Uncoupled  Yaw  Channel 
Autopilot;  Simplified  Estimator  Design; 
Discrete  Closed  Loop  System;  Elliptical 
Airframe 
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Figure  3-25  Yaw  Tail  Incidence  vs  Time;  Uncoupled  Yaw 
Channel  Autopilot;  Simplified  Estimator 
Design;  Discrete  Closed  Loop  System; 
Elliptical  Airframe 


The  pole-zero  plot  of  Figure  3-26  indicates  that  the 


discrete  closed  loop  system  is  stable,  since  the  z-plane  closed 
loop  poles  were  found  to  be: 


= 

0.0967149 

(III .C. 2-9) 

0.830876+j0. 31913 

(III. C. 2-10) 

% 

= 

0.830876-j0. 31913 

(III. C. 2-1 1 ) 

^4 

0.970259+j0. 0325175 

(III. c. 2-13) 

= 

0.970259-j0. 0325179 

(III. C. 2-14  ) 

= 

0.894774 

(III. C. 2-15) 

^7 

= 

0.937225 

(III. C. 2-15) 

The  time  response  plots  of  the  roll  angle,  angular 


rate  and  tail  incidence  are  presented  in  Figures  3.27  through 
3.29.  A  close  observation  of  the  above  pole-zero  and  time 
response  plots  for  the  discrete  roll  estimator  design  indicates 
that  they  are  very  close  to  those  of  the  discrete  classical  design 
found  in  the  previous  chapter. 

c.  Simplified  Design 

The  estimator  roll  designed  autopilot  can  be 
simplified  by  reducing  the  returning  gain  loops  as  follows: 

0 
0 
0 

0.0178  (III. C. 2-17) 

-0.0179 
-1.499 
0 

The  pole  and  time  response  plots  of  the  resulting 
simplified  estimator  roll  autopilot,  shown  in  Figures  3.30 


through  3.33,  do  not  present  significant  differences 

from  the 

corresponding  plots  of  the  discrete  classical  design. 

The 

discrete  closed  loop  system  is  again  stable,  since  the  z-plane 

closed  loop  poles  are: 

=  0.0992745 

(III .c. 2-18) 

^2.  -  0.836632+j0. 305719 

(III. C. 2-19) 

=  0.836632-j0. 305719 

(III. C. 2-20) 

^4  =  0.967647+JO. 0334948 

(III. C. 2-21 ) 

2^-  =  0.967647-j0. 0334948 

(III. C. 2-22) 

Zc  -  0.895937 

(III. C. 2-23) 

Z-j  -  0.938256 

(III  .C. 2-24  ) 

TIME  -  SEC 


Figure 


3.27  Roll  Angle  vs  Time;  Uncoupled  Roll  Channel  Autopilot 
Estimator  Design;  Discrete  Closed  Loop  System; 
Circular  Airframe 


TIME  -  SEC 


Roll  Angular  vs  Time;  Uncoupled  Roll  Channel 
Autopilot;  Estimator  Design;  Discrete  Closed  Loop 
System;  Circular  Airframe 
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IV. 


MODERN  CONTROL  DESIGN  AND  ROBUSTNESS  ANALYSIS  OF 
COUPLED  PITCH  AND  ROLL  CHANNEL  AUTOPILOT 
CIRCULAR  AIRFRAME 


A.  GENERAL 

The  present  chapter  deals  with  the  modern  control  design  and 
robustness  analysis  of  the  discrete  coupled  pitch  and  roll 
channel  autopilot  for  the  circular  airframe  configuration.  The 
continuous  open  loop  coupled  autopilot  whose  plant  system  and 
input  matrices  are  presented  in  Appendix  G  is  obtained  by 
coupling  the  linear  uncoupled  pitch  [Ref.  10]  and  roll  (Table  IV) 
channels.  Then,  utilizing  analog-to-digital  conversion  by  the 
aid  of  ORACLS  program  and  for  a  sample  period  of  0.0125  seconds, 
the  seventeenth-order  discrete  coupled  system  with  matrices  shown 
in  Appendix  H  is  formulated.  Next,  introducing  the  control-law 
and  estimator  designs  were  obtained  and  analyzed  in  terms  of  thei 
transient  responses  and  the  application  of  the  POPLAR  design 
program  [Ref.  7].  The  POPLAR  program  is  applied  in  order  to 
employ  singular  value  analysis  and  the  use  of  an  optimization 
routine  to  aid  in  pole  placement  control  design  of  the  above 
discussed  linear  multivariable  systems.  The  robustness  of  the 
system  is  also  considered  by  establishing  singular  value  levels 
which  correspond  to  multiloop  gain  and  phase  margins  determined 
from  the  universal  gain  phase  system  diagram  developed  by  Newsom 
and  Mukhapadhyay  at  NASA  Langley. 


B.  DISCRETE  COUPLED  STATE-FEEDBACK  DESIGN 
1 .  Design  Approach  and  Analysis 


The  discrete  coupled  state-feedback  designed  autopilot  is 
formulated  by  introducing  into  the  original  control  system  of 
Appendix  H  the  following  combined  pitch  [Ref.  10]  and  roll 
(III. B. 1-5)  control-law  gain  vector. 
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(IV. B. 1-1  ) 

The  pole-zero  plot  of  Figure  4.1  indicates  that  the 


discrete  coupled  system  is  marginally  stable,  since  the  z-plane 
closed  loop  poles  are: 


= 

0. 1  3446  +  jO. 0320974 

(IV. B. 

1-2) 

= 

0.13446-jO. 0320974 

(IV. B. 

1-3) 

^3  = 

0.89699+j0. 0906457 

(IV. B. 

1-4) 

0.89699-j0, 0906457 

(IV. B. 

1-5) 

■ 

0.955624+j0. 0299994 

(IV. B. 

1-6  ) 

= 

0.953624-j0. 0299994 

(IV. B. 

1-7) 

= 

/ 

0.998202 

(IV. B. 

1-8) 

= 

C'  ^ 

1 

(IV. B. 

1-9) 

-d,  = 

0 .995608  +  j0.083071  3 

(IV. B. 

1  -1  Q  ) 

0.995608-j0. 0830713 

(IV. B. 

1-11) 
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=  0.0993086 

(IV. 5. 1-12; 

^,2  =  0.836555  +  j0. 306107 

(IV.  3. 1-13) 

£,3  =  0. 836553-00.306107 

(IV. B. 1-14) 

=  0.967788+00.0334452 

(IV. B. 1-15) 

£■,<  =  0.967788-00.0334452 

(IV. B. 1-16) 

=  0.89573+00.938133 

(IV. B. 1-17) 

^,7  =  0.89573-00.938133 

(IV. B. 1-1 3 ) 

The  time  response  plots  of  the  roll  angle,  angular  rate 


and  tail  incidence  are  presented  in  Figures  4.2  through  4.4.  A 
close  observation  of  these  plots  indicate  that  they  are  identical 
with  those  of  the  discrete  uncoupled  state- feedback  design  of  the 
previous  chapter. 

2 .  Robustness  Analysis 

Executing  the  POPLAR  design  program  of  [Ref.  7]  with 
inputs  the  data  presented  in  Appendix  I,  the  minimum  additive 
input  (MIN  ADD  I.N  SV)  and  output  (SVADMO)  singular  values  were 
computed  from  a  frequency  range  from  0  to  200  rad/sec. 

Figures  4.5  and  4.6  which  are  plots  of  SVADMO  and  MIN  ADD 
IN  SV  versus  frequency  indicate  that  the  discrete  coupled  state- 
feedback  design  is  robust.  It  is  noted  that  for  very  low 
frequencies  the  values  of  SVADMO  are  above  0.8680. 

Finally,  in  terms  of  optimization  results,  the  ordered 
computed  eigenvalues  are: 

=  0.07717  (IV. B. 2-1 ) 

=  0.09855 

^5  = 


0 . 25923 


(IV. B. 2-2) 
(IV. B. 2-3) 
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Figure  4.2  Roll  Angle  vs  Time;  Coupled  Pitch  and  Roll  Channel 
Autopilot;  State-Feedback  Design;  Discrete  Closed 
Loop  System;  Circular  Airframe 
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Figure  4. 


3  Roll  Angular  Rate  vs  Time;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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Figure  4. 


4  Roll  Tail  Incidence  vs  Time;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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gure  4.5  SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design; 
Discrete  Closed  Loop  System;  Circular  Airframe 
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Figure  4.6  MIN  ADD  IN  SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  State-Feedback  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 


=  0.069508 

Zc  =  0.83707-j0.30l4 
=  0.83707  +  j0.30l4 
2ij  =  0.89844 
f-y  =  0.90064 
=  0.43996 
=  0.46207 

=  0.96562-j0. 03401 
Zii  =  0.96502  +  J0. 03401 
Z,;  =  0.99559-j0. 08307 
^,4-  =  0.99559  +  j0. 08307 
=  0.99820 


^.7  = 


1.06546 


(IV.B. 2-4  ) 
(IV.B.2-5  ) 
(IV.B. 2-6) 
(IV.B. 2-7) 
(IV.B. 2-8) 
(IV.B. 2-9) 
(IV.B. 2-10) 
(IV.B. 2-1 1 ) 
(IV.B. 2-12) 
(IV.B. 2-13) 
(IV.B. 2-14) 
(IV.B. 2-15) 
(IV.B. 2-16) 
(IV.B. 2-17) 


C.  DISCRETE  COUPLED  ESTIMATOR  DESIGN 


1  .  Design  Approach  and  Analysis 


Following  the  same  procedure  as  in  the  case  of  the  state- 
feedback  but  for  the  coupled  estimator  gain  vector  the  pole-zero 
and  time  response  plots  of  the  coupled  estimator  design  are 
obtained.  The  time  responses  are  again  identical  with  those  of 
the  discrete  uncoupled  estimator  design. 

2 .  Robustness  Analysis 

Executing  the  POPLAR  design  program  of  [Ref.  7],  Figures 
4.7  and  4.8  are  obtained  which  prove  the  robustness  of  the 
system.  It  is  noted  that  for  very  low  frequencies  the  values  of 
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Figure  4.8  MIN  ADD  IN  SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channel  Autopilot;  Estimator  Design;  Discrete 
Closed  Loop  System;  Circular  Airframe 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  goal  of  the  present  thesis  was  the  design  and  analysis  of 
discrete  lateral  autopilots  for  application  to  BTT  missiles.  The 
following  are  the  principal  conclusions  based  on  this  work. 

1.  The  continuous  and  discrete  classical  designed  autopilots 
were  proved  to  have  identical  performances  for  the  two 
lateral  channels. 

2.  The  state- feedback  and  estimator  autopilots  were  introduced 
as  additional  dynamic  designs  in  order  to  implement  control 
to  the  original  system.  Both  designs,  analyzed  in  terms  of 
their  transient  responses,  were  found  to  meet  the  desired 
requirements . 

3.  The  simplified  state-feedback  and  estimator  designs  reduced 
some  of  the  returning  gain  loops,  making  the  system 
simpler,  without  any  significant  effects  on  the  system's 
performance . 

4.  The  performance  of  the  coupled  pitch  and  roll  channel 
autopilot  was  found  to  be  satisfactory  and  the  overall 
system  proved  to  be  robust. 


B.  RECOMMENDATIONS 

In  order  to  improve  the  simplicity  of  the  overall  system, 
more  returning  gain  loops  of  the  state-feedback  and  estimator 
design  must  be  eliminated.  A  further  investigation  then  must  be 
conducted  in  order  to  examine  if  the  performance  of  the  resulting 
design  remains  unchanged. 


APPENDIX  A 


DESIGN  REQUIREMENTS  FOR  UNCOUPLED  AUTOPILOTS 

The  requirements  for  the  classical  design  method  of  the 
uncoupled  channel  autopilots  [Ref.  6]  are  the  following: 

1 .  High  Frequency  Attentuation  in  Actuator  Command  Branch 

a.  Uncoupled  Yaw  Channel 

It  must  be  ^  15  db  at  100  rad/sec  and  zero  angle-of- 

attack  and  sideslip.  This  requirement  will  provide  sufficient 
high  frequency  attenuation  for  ^  30  Hz  actuator  and  for  body 
bending  modes  when  high  frequency  filters  are  added,  but  it 
limits  the  ability  of  the  yaw  autopilot  to  minimize  sideslip 
angle . 

b.  Uncoupled  Roll  Channel 

It  must  be  >y  15  db  at  100  rad/sec  and  zero  angle-of- 
attack.  This  will  provide  sufficient  high  frequency  attenuation 
for  30  Hz  actuator  and  for  elastic  modes  when  high  frequency 
filters  are  added,  but  this  requirement  limits  the  speed  of  roll 
angle  response. 

2 .  Relative  Stability  for  Both  Lateral  Channels 

Gain  margins  >6  db,  phase  margins  ^30°  with  a  goal  of 
12  db  and  50°. 

3 .  Acceleration  Time  Response 

a.  Uncoupled  Yaw  Channel 

(1)  63^  time  constant  of  approximately  0.4  seconds. 
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(2)  Overshoot  10?. 

(3)  Steady-state  error  need  not  be  zero 
Uncoupled  Roll  Channel 

(1)  63$  time  constant  of  0.5  seconds. 

(2)  Overshoot  10$. 

(3)  Zero  steady-state  roll  angle  error. 


APPENDIX  B 


AERODYNAMIC  DATA 

The  overall  classical  design  developed  by  Arrow  [Ref.  6]  was 
performed  for  the  selected  flight  condition  of  Mach  number  3.95 
at  60000  feet  altitude.  The  corresponding  aerodynamic  data 
presented  below  were  taken  or  derived  from  the  ICAO  standard 
atmosphere  tables. 


TABLE  VII 

AERODYNAMIC  DATA  (M=3.95,  H=60kft) 


Temperature,  T  (®R) 

389.988 

Sonic  Velocity,  a  (ft/sec) 

968.47 

Pressure,  p  (Ib/ft^) 

149.78 

Density,  (Ib-sec^/ft) 

0.0002238 

Velocity,  V  (ft/sec) 

3825.4565 

Dynamic  Pressure,  q  (Ib/ft^) 

1637. 145 

APPENDIX  C 


LINEARIZED  AERODYNAMIC  DERIVATIVES 


The  linearized  aerodynamic  derivatives  at  the  selected  flight 


condition  and  about  a  zero  trim  angle-of-attack  are  provided 


below  for  both  airframe  configurations  [Ref.  6]. 


TABLE  VIII 


LINEARIZED  AERODYNAMIC  DERIVATIVES 
(M=3.95,  H=60kft;  3=0°) 


CIRCULAR 


-0.065 


-0.025 


0.021 


-0.050 


0.031 


ELLIPTICAL 


-0.043 


0.024 


0.016 


-0.042 


0.023 


*»'  ^ 
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APPENDIX  D 

MISSILE  SIZING  AND  MASS  PROPERTIES 

In  order  to  provide  a  realistic  missile  based  on  the 
aerodynmamically  tested  configuration  concepts  [Ref.  5],  the 
models  were  assumed  to  be  1 /6-scale  and  the  mass  properties  were 
developed  corresponding  to  mass  distribution  which  might  be 
expected  for  missiles  of  this  size.  All  the  geometric  and  mass 
properties  are  presented  in  the  following  table. 

TABLE  IX 

GEOMETRIC  AND  MASS  PROPERTIES  OF  MISSILE  CONFIGURATIONS 


Length,  1  (in) 

Max.  Diameter  (in) 

Max.  Major  Axis  (in) 

Max.  Minor  Axis  (in) 
c.g.  distance  from  L.E.  (in) 
Reference  Length,  d  (ft) 
Reference  Area,  S  (ft) 
Weight,  W  (lb) 

(slug-ft^) 

1^2  (Slug-ft^) 


Circular 

168 

24 


100.8(0.6  ) 
2 
TT 

2525 

40 

810 


Elliptical 

168 

41.57 

13.86 

100.8(0.6  ) 

2 

TT 

2475 

40 

853 
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APPENDIX  E 


PROGRAM  LOGIC  FOR  APPLICATION  OF  ACKERMANN’S  FORMULA 

The  program  logic  for  computing  the  control-law  gain  vector  F 
or  the  transpose  of  the  estimator  gain  vector  taken  from 
[Ref.  9],  is  given  in  the  following  table. 

TABLE 

PROGRAM  LOGIC  FOR  APPLICATION  OF  ACKERMANN’S  FORMULA 


I.V 

14. 

1. ^. 
Ih. 
17. 
IS 
IV 
2(1. 
21. 

23. 

24. 

2. “'. 
26. 
27, 

25. 
2V 
■'(). 
31 


Re.Kl  in  <t>.  I'.  T.  ;inj  the  numher  ol  ".(ates. 

Comment:  lir>l  we  will  read  in  the  dcMied  pole  locations  m  the  \-pliine.  convert  them 
to  :-pUine  poivnomial  coct'hcients.  and  construct  ai‘l>). 

I  —  identity  matrix.  \,  »  S, 

M.PIIA  —  I 

k  —  I 

If  k  •  V, .  ito  to  Step  IS. 

Head  in  pole  location  C  us  «  ~  jh. 

If  h  -  II.  CO  to  step  14. 

.1,  *-  -  2  cxpii/7  I  cos  hT 
•tj  —  e\pi2i(7  I 

\LI‘M\  —  Al.l’IlA  .  i<|.  ■  <1*  -  t  .<|>  -  I, II 
k  —  k  -  2 
Go  to  step  6 
.4 ,  evpi.i  /  I 

AI  I>M  \  —  M.rilA  '  l<l«  -  A,  '  n 

k  ^  k  -  \ 

Go  to  step  6. 

Comment:  now  we  construct  the  controlhihiliiv  matrix. 

C  —  I 

E  -  r 

k  —  I 

If  k  •  .N  . .  ito  to  step  2S. 

Ci'mment:  replace  column  k  of  C'  h\  E. 

Cl  :  A  I  ^  F. 

A  —  A  *  I 
K  —  'l«  E 
Go  lo  step  22. 

Comment:  now  solve  for  the  contiol  l.iw  .  lust  form  ('as  the  last  row  of  I. 

1 

.SoKe  HC  =  E  for  B. 

K  -  B  ■  AI.IMIA 

i-M) 
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APPENDIX  F 

ACKERMANN  FORTRAN  PROGRAM 


LE:  ACKERMAN  mAXrlV  Al 


integer  NS.K.U,  IEa.1  fJ.NI 

real  Ph(20.2U I ,GA<  20,2 JI.au  20,  20) , ALPHA (20, 20) , nKl 40) 

REAL  A2  .E  (20.2  J)  ,0(  20.20)  ,G(20)  .  EE  1,:J, 2  0)  .JK(  20,2  J) 

REAL  A.o.T,  R,ATEME(20,20) 

REAL  Hi  , TU120,2Q1 , TK 20,201 . TS ( 20.20) .0(20,20 ), TM(20 ,20) 

INPUT  OF  NS=  NUMBER  OF  STATES 
T=  Sample  Hme 

PRINT  ,  •  ENTER  NUMBER  OF  STATES.* 

REAJ,  NS 

PRINT,  •  ENTER  SAMPLE  TIME  • 
real),  T 

INPUT  JF  PHI=A  MATRIX 
00  lO'J  1  =  1,. NS 
DC  101  J=l,.',S 
RRITe(6,I03)  I.J 

103  FORMAT!  4X,  ■  £N  ie«  P  H(  *  .  12  .  *  .  *  ,  12,  *  )  *  ) 

REA0.Ph( I , J ) 

101  continue 
100  continue 

INPUT  OF  GA=0  MAT.RIX 
00  200  1=1 , NS 

J  =  l 

kHIT£(6,105)  I.J 

105  FORMAT!  4X,  •  ENTEa  GA (  •  ,  12  .  * ,  *  ,  12 ,  *  )  *  ) 

REAO.GA (I ,J ) 

200  CONTINUE 

aUlLOING  THE  lOENTITY  MATRIX 
AlPMA=C= I 

00  1  1=1. NS 
DC  5  J=  I, NS 
IF  (  I  .£  C.J  )  GO  TO  300 
Al(  I,J)=0. 

ALPHA!  I  ,J)  =  0. 

C!1  ,U)  =  0. 

GC  TO  5 

300  AU  I.J)=1. 

ALPHA!  I  ,J)  =  1  . 

C!1,J)=1. 

5  CONTINUE 

1  continue 

K=1 

16  IF  (K.GT.NSJ  GO  TO  U 

INPUT  OF  CEiUEC  POLE^  LOCATION 
hRITE  16.43) 

43  FURMAI!  •  ENTER  REAL  PARI  OF  OESIRE'J  POLE  LOCATION.*) 

READ, A 

MRlTEi6.44) 

44  FOR, MAI! ■  enter  IMAGINARY  PART  OF  DESIRED  POLE  LOCATION.') 

REAO.d 

BUILDING  alpha  (PHI) 

IF  (a  .EJ.J)  go  to  14 
Al=-2.»exP(  h*T)  *COS(d*TI 
42=exp( 2. =4*1 ) 
evaluate  PF1*PHI 

CALL  VMULFF  !  PH  ,PH  ,  NS  ,  NS  .  NS  ,  20  ,20  ,  TM  ,2  0,  ICR  ) 

DC  2  I =1,NS 
OJ  6  J=  l,NS 

Toil  .  J)  =  PH(  I,  J)»41 
Till,  J|  =  AI  !  I,  J)*A2 
T  j!  I  ,  J)  =  rM!  I,  J)*IU!  I  ,J  )  *11  (  I  ,J) 

6  CCNIINJE 

2  CCNTINUE 

CALL  VMULFF  I  AL  PHA  ,  T  S .  NS  .  NS  ,  NS  ,20 , 2  0  ,A  TEMP,  20  ,  I  E  R) 

UO  600  1=1. NS 
00  601  J=1.NS 
601  ALPHA! I ,J|=4TEMP! I ,U| 

6C0  CONTINUE 


AOKOOOlO 
ACKO0O2O 
ACkGCOBO 
AC  kOOOhO 
ACkuOOSO 
ACKCU060 

ACK  joo  ro 
AS  KULiO  50 
ACK0C090 
ACKOOlOO 
ACKOCllO 
A0KCO12O 
ACKOul 30 
ASKOOl 40 
AS, <001  50 
ASKC.lbO 
ASXO'SITJ 
ACkOJI dO 
ASxOolHO 
ASh00200 
AS<ao2 10 
ASKJS220 
AS<0C230 
ACK00240 
AC<0U250 
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APPliNUIX  11:  PLANT  SYSTLM  AND  INPUT  MATRICES  OF  DISCRETE  COUPLED  PITCH  AND 
ROLL  CHANNEL  AUTOPILOT 
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APPENDIX  I 


COUPLED  STATE-FEEDBACK  DESIGN  INPUT  DATA  FOR  POPLAR  PROGRAM 


FILE:  FEEDBACK  DATA 


A1 


1000  0 

200.0  0 

1 

0.4  0 

0  0  0 

17171702171702171 7 
1.5  /y5758E-01 
1.43769Q5E-03 
0 
0 

7.3702121E-Q3 
6.48032  77E-C*> 

0 

0 

-3. 5755Q72E-2 
-3. 14630456-5 
0 
0 

6.36419146-04 

1.03a917aO£-7 

0 

0 

-1.2632605e-» 

-l.l935039E-a 

0 

0 

5.00535036-05 

9.9653837E-01 

0 

0 

1.  1 1980  72£-02 
-5.53283156-1 
0 
0 

-1.8197865E-3 

-7.69372686-7 

0 

0 

-2.  4201  113E-3 
-2.20856126-6 
0 
0 

3.51203616-03 
3.46475186-03 
0 
0 
0 


1.0 

0 

0.0  0 

3  3  1000 


3.63693266-02 
-o. 85501956- 3 
0 
0 

1. 00018476-00 
-4. 31433676-5 
0 
0 

—  a.  968  35566-4 
2 .094o996c-04 

0 

O 

4.30723476-03 

—  1.  02  769066-0 

0 

O 

-3.31021336— * 
1.34677576-37 
0 
0 

3.49067026-04 
1  .24738786-02 
0 
0 

8.45897156-02 

9.94659846-01 

0 

0 

-1.64759366-2 
o  .11292606-06 
0 
0 

-2.62507036-2 

l.to5903yE-05 

0 

0 

7.717026-66-02 

—  1. 25026406-2 

0 
0 
0 


5.58485  926 -03 
3.  121503806*0 
0 

1.21457136-02 

2.35544156-02 

0 

9.39275546  -01 
-1.14342466-1 

0 

7.37304526 -04 
-2.57727336-1 
0 

-1.35084526— » 
1.10524796-01 
0 

5.35213356-05 

-4.383796oE-2 

0 

1.29826646-02 
-4. 531317.;  6-0 
0 

-2.53428716-3 

1.9O53256C-02 

0 

-"t.04o8  1  1 1  6-  3 
-l.Ol J45586-1 
0 

1.  18702  99E -02 
4.65o3462£— Ou 
0 


-3. 54l5236£-2 
-2  .  74  3‘*6  95  6  *0 
0 

-1.31103296-4 
-2.  16579256-2 
0 

3.79232406-04 
1.05  1  a*  9'*c-0 1 
0 

9.91345926-01 
-*♦•  37^06766  —  1 
0 

-5. 52425915-1 
1.00735536-01 
0 

-3.42352576-4 

-4.31031496-2 

0 

-3.25733366-2 

-7.6933073E-0 

0 

1. 59 144376-02 
8.25513205-01 
0 

2.50737846-02 
8 .52133306-01 
0 

-7. 453506OE-2 
-3.  84^95976-0 
0 


-5.63506346-3 

8.496u3416-Cl 

0 

-2.37079036-5 

9.078»743£-03 

0 

1.3937755t-04 
-•*  .4066  55  1  c  -2 
3 

1.17240 94C-02 
3.02396946-04 
0 

9.94314316-01 
-5 .6  5300o2  6-5 
0 

-5.42321945-5 

2.239022306-5 

0 

-1. 3194..30E-2 
6.37744246-03 
0 

2.59274996-03 

-1.61176916-3 

0 

4.  16749  756-03 
-3  » Oo06  4  3  76  —  3 
0 

-1.21446506-2 

I.0o75ol0c-00 

0 


0 

0 

0 

0 

3 

9.97738786-01 

3.22851556-00 

-2.  7C1  36906-3' 
3.o730433e-02' 

1.02458966-03 

2.26857366-01 

-2.23924736-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

1.249  26126-02 
1.61206056-02 

9.S9S'Vl2U-0l 
3«2i^5b56— J4 

3.36863396-06 

7.92138316-04 

-7.32510196-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.32996046-3 

-9.65973186-4 

-2.  09357296-  1 
-2.O900370E-5 

9.04837276-01 

-3. 68819656-5 

3.64511 286-05 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-9.  07977236-3 
-1.70C7972E-1 

-  1.  1C2408  76-2 
-2.  00285126-3 

3.99972926-03 

9.27630326-01 

1.16307326-02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1. 73409426-4 
2.45055106-01 

-2.07178996-4 

2.725751U6-03 

7.85645386-05 

1.73740496-02 

9.10594  316-01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.90810306-3 

7.46927476-01 

-2.  30a  18306-3 
-  1.  142948  56-3 

8*4o098U5c— 04 

1.52307156-01 

-1.50374136-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

156 


i  lOOMd  JH 


file:  FEEOdACK  JATA 


-6.287571Bt-2 
4. 3323648e»Cl 
-7.  5609329E-3 
-1.64491B6E-2 
7. 932162  3c-02 
-9.9881608E-4 
1.8285612E-04 
-7. 2448 726C-5 
-1. 7575515E-2 
3. -.31  55566- C3 
5.4S07555E-C3 
-1.O073019E-2 
0 
0 
0 
0 
0 
0 
0 

i.oacaooae-co 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

u 

0 

J 

0 

!•  i3  COCO  OOE  —  oO 
0 
0 
0 
3 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
a 
0 
0 
0 
0 
0 

l.OOOOOOOE-CO 

0 

0 

0 

u 

0 

0 

o 

0 


7.  5602497E-2 
.13406856-02 
0 
0 
0 
0 
0 
0 
0 
0 
Q 
0 

.7C1364DE-03 
.792o971r-0o 
.39357296-31 
.10240a7£-Q2 
.07173996-34 
.3081 830E-03 
.5oQ2497E-02 
0 
0 
0 
0 

.OCOOOOOE-JO 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.COCUOOOE-DO 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.00000006-30 

0 

0 

0 

0 


2.3l7ul9l-3-02  5.O1224356-00  -5.53140556-0 


1 .00000006-00 
0 
0 


1.00000006-00 

0 

0 


1.00303006-0 


1 .COOOOOOE-OO 
0 


1 . 00000 006-00 
0 


1 .00003  U06-C 


1 .00000006-00 


ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo  ooo 


ILE:  FEEOdACH  0AT4 


A  1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.OO0QOOUE-C3 

0 

0 

0 

J 

-o.oiaaoooE-oo 
O.OUOOOOOEtCO 
0.  OOOdUOOEi-QO 
O.QQGOOOOEi-CO 
0. JOOOUOOE+OO 
C.OaOOOOOb-GO 
C.OCOOOOOfc*-CO 
-1.4S9QQUUE-dO 
0  Q  0  0  0  0 
0  0  0  0  0  0 
-lOOQO 
0.13199300 
0. 13199300 
0.39698700 
0.89698700 
0.95362500 
0.95362500 
0.99  82  02  00 
l.OOOQQOOO 
0.99560800 
0.99360800 
0.0992 8050 
0.33521600 
0 .33321600 
0.96903700 
0.9o908700 
0.39379700 
0.93707500 


O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1  .OCCQOOQE-JO 
0.0279000E-C0 
O.OOOlOOOctOO 
U.OCOOOJUEfOO 
0. JOOOOOOE+O 
a.uccj000E«’00 
O.QCCOOOuc+Ou 
0 . 0  wO'.OOOc  —00 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0042300E-03 
1.7732000E+00 
O.OOOOOOot  *-00 

Q.COOOOOOE  *00 
C. 0003000 E»0 
0-C300000t»00 


0 

I .OOOOOOOE-uO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 02900 00d*00 
2.016200JE+00 
0.  OOoOOOOE-^OO 

0.00000006*3 
0. OOOOOOjE*0 
0.01 78000E-0 


Q.0C23Q00E-OO 
000  OOCOOQOO 
COO  OOCOOOOO 
10000 

0.03  1937  70 
-0.03193770 
0.09067030 
-G.0Svj670dQ 
C.02  9985UQ 
-0.02998500 
Q.OQOOOOOO 
C.UOOCCOOO 
0.033071 30 
-0.033071 30 
O.OOOOOOoO 
C.  30  >73630 
-0.309  73600 
C. 03235880 
-0. 03235880 
C.OOOOCOOO 
C. 00000330 


0 

0 

0 

0 

1.03003006-00 

0 

0 

0 

0 

0 

0 

o.coAiooce+oo 
o.uijio Joe*oo 

0.0000000fc*00 

0.30033306*03 
0 . 03u3  OOOE  *0 
-0.01 79000E-0 


APPENDIX  J 


COUPLED  ESTIMATOR  DESIGN  INPUT  DATA  FOR  POPLAR  PROGRAM 


FUE:  cSTlMAfO  JmTA 


AL 


lOGO  0 

100.0  O.G 

1  0 
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